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The results of recent phenomenological studies of unpolarized nuclear deep-inelastic scattering are discussed and 
applied to calculate neutrino charged-current structure functions and cross sections for a number of nuclei. 



1. Introduction and motivations 

Significant nuclear effects were discovered 
in charged-lepton (CL) deep-inelastic scattering 
(DIS) experiments (for a review see |1I2| ). The 
experimental observations indicate that the nuc- 
lear environment plays an important role even at 
energies and momenta much larger than those in- 
volved in typical nuclear ground state processes. 
The study of nuclei is therefore directly related 
to the interpretation of high-energy experiments 
with nuclei from hadron colliders to fixed target 
experiments. 

The understanding of nuclear effects is partic- 
ularly relevant for neutrino processes in which 
weak interaction with matter requires the use of 
heavy nuclear targets in order to collect signifi- 
cant number of interactions. Therefore a reliable 
treatment of nuclear effects is important for in- 
terpretation of neutrino experiments and in some 
cases crucial for reducing systematic uncertainty. 

In this paper we briefly review the results of re- 
cent studies of CL nuclear DIS of Ref. [3] and ap- 
ply this approach to calculate the (anti)neutrino 
differential cross sections of charged-current (CC) 
interactions with nuclei. 

2. Nuclear structure functions 

The theoretical background of the analysis of 
Ref. [3] involves the treatment of a few different 
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mechanisms of nuclear effects which arc charac- 
teristic for different kinematical regions. Summa- 
rizing, for the nuclear structure function of type 
a = 1, 2, 3 we have 

= K"" + K'^ + K'^ + Scoi.F^, (1) 

where Fa^^ [Fa^^) are the incoherent contri- 
butions from the bound protons (neutrons) cor- 
rected for Fermi motion, nuclear binding and off- 
shell effects (impulse approximation). The term 
Fa^^ is a correction associated with scattering off 
the nuclear pion (meson) field. The term 6co\iF^ 
is a correction due to coherent interaction of in- 
termediate virtual boson with nuclear target. 

The terms Fa^^ and Fa^^ are calculated as the 
proton and the neutron structure function aver- 
aged with the nuclear spectral function and cor- 
responding kinematical factors. In particular, for 
F2^^ we have^ 

Fl'^{x, Q2) ^ ((1 + k,/M)F^{x\ Q\ e)) , (2) 

where k = {M + e,k) and x' = x/{l + {e + kz)/M) 
are the four-momentum of bound proton and its 
Bjorken variable. The averaging is taken over the 
proton spectral fmiction which describes the dis- 
tribution of bound nucleons over momentum and 
energy. Similar expression holds for the neutron 
term. 

^Eq.J2J is written for the kinematics of the Bjorken Umit 
assuming that momentum transfer is along z axis q = 
(go, Ox, — |q|). For more general expressions valid at finite 
Q and further detail see Ref. 1^ . 
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Nuclear Fermi motion and binding corrections 
provide a correct trend of observed behavior of 
the ratios 112 = F^/F^ (so-called EMC ratio) 
however the quantitative description of data is 
missing in impulse approximation. In a quan- 
titative treatment it is important to go beyond 
this approximation and take into account that 
the structure functions of bound nuclcons can 
be different from those of the free nucleons. In 
Ref.[3] this effect is modelled by off-shell correc- 
tions, i.e. the dependence of structure functions 
on fc^ in Eq. (j^l . Since characteristic energies and 
momenta of bound nucleons are small compared 
to the nucleon mass we treat off-shell effects as 
a linear correction in fc^ — to the structure 
function of the on-shcll nucleon 



F2{x,Q\e)=F2{x,Q')[l + Sh^ 
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The function df2{x,Q^) describes the relative 
off-shell effect. In analysis of Ref.|3] this func- 
tion is treated phenomcnologically. In particu- 
lar, we assume 6/2 to be independent of and 
parametrize its x dependence as 



<5/2 = Cn{x - xi){x - xo){h - x), 



(4) 



where < xi < xp < 1 and h > 1. 

It should be noted that because of binding the 
nucleons do not carry all of the light-cone momen- 
tum of the nucleus and the momentum balance 
equation is violated in the impulse approxima- 
tion thus indicating the presence of non-nucleon 
degrees of freedom in the problem. In Ref. [H] we 
consider a correction to nuclear structure func- 
tion due to scattering off meson fields in nuclei in 
the convolution approximation 
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= / dyU/Aiy)F^{x/y,Q'), 
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where fT^/Aiv) is the distribution of nuclear pion 
excess in a nucleus and F2 is the pion structure 
function. The distribution function /tt/^ (y) is cal- 
culated using the constraints from equations of 
motion for interacting pion-nucleon system. By 
using the light-cone momentum balance equation 
we effectively constrain the contribution from all 
mesonic fields responsible for nuclear binding. 



In the small-x region the coherent effects in DIS 
are relevant (for a review see Ref.|2]). These ef- 
fects are associated with the fluctuations of in- 
termediate virtual boson into quark-gluon (or 
hadronic) states. At small x an average time of 
life of such fluctuation is significantly larger than 
the average distance between bound nucleons. 
For this reason the virtual hadronic states un- 
dergo multiple nuclear interactions while travers- 
ing a nucleus that causes nuclear shadowing ef- 
fect. The rate of this effect depends on the scat- 
tering amplitude of the virtual hadronic states off 
the nucleon. In our approach we model the inter- 
action of virtual hadronic states with the nucleon 
as scattering of a single state with effective cross 
section a. In this approximation the relative nuc- 
lear correction to the structure function is deter- 
mined by the corresponding correction to effective 
cross section d. In our studies PJ the effective 
cross section a for the scattering off the nucleon 
is treated phenomcnologically and parametrized 
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where the parameters ctq and tri describe low- 
Q and high-Q limits while the scale Qo controls 
transition region. The nuclear corrections to ef- 
fective cross section were calculated using the 
Glauber-Gribov multiple scattering theory. 

We used the outlined approach in the analy- 
sis of data on the ratios 'R2{A/B) = F^/Ff 
of structure functions of two nuclei. The gen- 
eral goal was to develop a quantitative model of 
nuclear structure functions which, from one side, 
would include the major mechanisms of nuclear 
scattering and, from the other side, would de- 
scribe the existing data with acceptable accuracy. 
We analyze data on 7^2 for a variaty of targets 
from D to Pb for a wide kinematical region (for 
more detail see Table 1 of Ref. P] ) . From prelim- 
inary analysis of data we observed strong corre- 
lations between some of the model parameters. 
In order to reduce the number of free parameters 
we used additional constraints. In particular, the 
parameters h and xq turned out to be fully cor- 
related and related &s h = \ -\- xq. The parameter 
xi is strongly correlated with Cjv and Qq. We 
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performed several fits with diflFerent fixed values 
of Xi and found that xi = 0.05 corresponds to 
the lowest and provides a good cancellation 
between shadowing and off-shell correction to the 
normalization of nuclear valence quark distribu- 
tion. From preliminary fits the best fit value of ai 
was consistent with zero and we fixed cti = in 
the final fit. The parameter uo was fixed to 27 mb 
(averaged meson-nuclcon total cross section in 
the vector meson dominance model \412\ ) in or- 
der to reproduce the photoproduction limit. The 
remaining parameters Cat, xq and Qo were ad- 
justed to reproduce data. The global fit of Ref.p] 
to all data results in Cn = 8.1 ± 0.3 ± 0.5, xq = 
0.448±0.005±0.007, Ql = 1.43±0.06±0.2 GeV^ 
with x^/d.o.f. = 459/556 (the last error is the es- 
timate of systematic or theoretical uncertainty). 
In order to test the model we performed a num- 
ber of fits to different sub-sets of nuclei in the 
region from "^Hc to ^°^Pb. The results are com- 
paitible within the uncertainties with the result 
of the global fit thus indicating an excellent con- 
sistency between the model and the data for all 
nuclei. 

3. Neutrino-nucleus inelastic cross sections 

In this paper we apply the developed model of 
nuclear structure functions to calculate the dif- 
ferential cross sections of (anti)neutrino inelas- 
tic scattering. In contrast to the CL scatter- 
ing, which is driven by electromagnetic interac- 
tion, neutrino interactions are characterized by 
the presence of both the vector current (VC) and 
the axial current (AC) contributions. The VC- 
AC interference gives rise to P-odd and C-odd 
terms in the cross section which are described 
by the structure function F3, which is absent 
in CL scattering. In contrast to VC the AC is 
not conserved. For this reason the AC contri- 
bution plays important role and even dominates 
the (anti)neutrino cross sections at low Q. In 
this region the terms due to non conservation of 
the AC are most important for the longitudinal 
structure function Fl which is determined by the 
divergence of the current (Adler theorem 0). At 
low momentum transfer the divergence of the AC 
is dominated by the pion field (PCAC) that al- 



lows us to calculate the leading term at low Q^: 
Fl'^'^^' = fT^aT^{s,Q'^)/'K, where /^r = 0.93 is 
the pion decay constant and (Ttt is the total pion- 
nucleon (nucleus) cross section with the center- 
of-mass energy squared s = M"^ + Q'^{l/x — 1). It 
should be also noted that cr^ describes interaction 
of an off-shell pion with momentum q and the vir- 
tuality —Q^. Although the contribution from the 
VC to Fl vanishes at low Q^, it rises with and 
has to be taken into account at ~ 1 GeV^ as 
well as for higher values of momentum transfer. 
We explicitly separate the PCAC term and write 
the full Fl as 

Fl{x, Q2) ^ i^r'VpcAc(Q') + Fl{x, Q'), (7) 

where Fl incorporates the contributions from 
VC and non-PCAC terms from AC which van- 
ish at 0. In order to interpolate between 
low and high we introduced a form factor 
/pcAc = (1 + Q^/MpcAc)"^ where MpcAc repre- 
sents the scale controlling the PCAC contribu- 
tion. Since the pion pole does not contribute to 
the structure functions (because of transversity 
of the lepton current, see Ref.[n] and references 
therein) the scale MpcAc is not determined by the 
pion mass but rather related to higher mass states 
like fli meson and pir continuum.'^ The term F^ 
dominates at high at which it is identified 
with the pQCD structure function with target 
mass (TMC) and higher twist (HT) corrections. 
In numerical applications we use the PDFs and 
HT of Ref. [3 which were derived from the anal- 
yses of CL DIS data. In order to evaluate Fl at 
low values of we apply polynomial extrapola- 
tion between pQCD calculation at high and 
the limit of Q — > with the matching point at 
Q2 = 1 GeV^ 0. 

The magnitude and the dependence of the 
PCAC term in F2 is illustrated in Fig. [H Note 
that the magnitude of the PCAC term decreases 
for heavy nuclei because of the nuclear shadowing 
effect for the pion cross section. 

It is also important to note that for neutrino 
scattering the PCAC term leads to the rising ratio 

■^In numerical calculations described below we used 
AfpcAc = Tiaj. In Ref.|H| it was argued that the rele- 
vant scale is given by the pir cut rather than the ai pole. 
However, both values are close numerically. 
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R = Fl /Ft at low in contrast to the CL case. 
This effect is illustrated in Fig. El where R{Q^) 
was calculated for two fixed values of x for the 
isoscalar nucleon (average over the proton and the 
neutron) and a number of nuclei |8I9| . 
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Figure 1. The PC AC term in neutrino F2 calcu- 
lated as a function of for the fixed x ~ 0.0001 
0. The results are shown for the isoscalar nu- 
cleon (average over the proton and the neutron) 
and for ^^C, ^^Ye and ^ospb nuclei. 

In Figures 01 El and El we report the results 
of our calculations of (anti)ncutrino differential 
cross sections in comparison with the measure- 
ments for several values of the (anti)neutrino en- 
ergy from 20GeV to 170 GeV 0. The calcula- 
tion of (anti)neutrino nuclear structure functions 
uses the same NNLO PDFs and HT terms ob- 
tained from the analyses of charged-lepton DIS 
data |7] and which were applied in the analysis 
described above. The treatment of nuclear effects 
is based on the results of Ref.[SI. In the calcula- 
tions we used the off-shell function (5/2, derived 
from CL data in Ref.[2|, for both F2 and F3 for 
neutrino and antineutrino. Nuclear effects in the 



PCAC term are controlled by the multiple scat- 
tering corrections to the pion cross section ct^. 
We also remark that nuclear pion excess correc- 
tion (Eq.©) vanishes for F3. As a result, we ob- 
serve a good agreement between the data and our 
calculations for all examined nuclei that provides 
a good test of the model of Ref.pi. It should be 
noted that the data points at low x bins, which 
typically have low (for example for CHORUS 
data ^ 0.25 GeV^), are also reproduced by 
calculation. In this region the cross sections are 
dominated by the PCAC term. Thus our analysis 
supports the presence of significant PCAC effects 
in neutrino inelastic scattering at 3> rn^. 
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X = 0.045 




x = 0.125 




Figure 2. Comparison of the ratio R = Fl/Ft 
calculated for charged leptons (dashcd-dotted) 
and neutrino (solid) for the isoscalar nucleon as 
a function of at a; = 0.045 (upper panel) and 
X = 0.125 (lower panel). Also shown are the re- 
sults of calculation of R for neutrino scattering 
off i^C, 56Fe and ^ospb (from top to bottom). 
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Figure 3. Comparison of our calculations (open 
symbols) with NOMAD data [TU] (full symbols) 
for neutrino differential cross-sections on ^^C at 
£■ = 20 GeV (upper panel) and £; = 60 GeV 
(lower panel). 



6 



S 

U 

CO 

9 1 

O 

i 

>* 
■D 

i 



0.2 
0.1 



>x = 0.020 



x = 0.080 




X = 0.045 



X = 0.125 



X = 0.225 




: X = 0.650 



0.2 0.4 0.6 0.8 1 

y 



0.2 0.4 0.6 0.8 1 

y 

• CHORUS V 

MODEL V 
A CHORUS anti-v 
^ MODEL anti'V 



; 0.78 
0.5 
0.25 

0.2 

0.1 
























r 


























— *— , 














■ , , , 


, , , 


, , , 


x = b.015 

, , 1 , , 








X = 0.045 

, , 1 , , 




























X = o.bso 








x = 0.125 


" , , , 


, , , 


, , , 












, 




*- 






4— S 














" , , , 


, , , 


, , , 


X = U.1 IS 

-A ^ 

, , 1 , , 








X = [j.zza 

1 , , 


t- 






1 














" , , , 


, , , 


, , , 


X = 0.275 

^T* i , , 








X = 0.350 

"1 , , i , , 


















» - n ssn 








X = 0. 


450 






































X = 0.650 








X = 0.750 


■ , , , 


, , , 


-A 

, , , 








-i-*- 


4 







0.2 0.4 0.6 0.8 II 

y 



0.2 0.4 0.6 0.8 1 

y 



..r=f=H^-i^H ■■■■ 

X = oi)20 



"-*^--i---f 

x = 0^)45 



S of 

2[ 

u 

? 1 

o 

i 01 
■o 

i 1 




0.2 





iX = 0.550 



t-r-t-.-i-. 



0.1 

oi 



'1 '"i 



IX = 0.650 

:t-t-( ..; 



0.2 0.4 0.6 0.8 1 

y 



0.2 0.4 0.6 0.8 1 

y 

• CHORUS V 

MODEL V 
» CHORUS anti-v 
- MODEL antl-v 



> 



b 1 

X 

1 1 



LU 



0.78 
0.5 
0.25 

0.2 

0.1 






































— s- 


-- » 


■ , , , 






x = 0.015 

. . i . . 








. - 4 ... 

X = 0.045 

. . i . . 








-|-"->-'4-, ■■■■ 

X = 0.080 








X = 0.125 


- , , , 






-A — ^_ 












-4 




T »- 


-4 1- 
















- t-i 












— 










X = U.I Its 
' » : * 








X = 0.225 

^-^-^ , 






H — r 


-|-+- 


-» 












" , , , 






X = 0.275 






~1 — 


x = 0.350 

-t-7 » 1 , , 


















» - n 




* 




-^^v 

y - n Afin 






















, , 






--♦ ■■■■ 






























X = 0.650 








x = 0.750 























0.2 0.4 0.6 0.8 II 0.2 0.4 0.6 0.8 1 

y y 



Figure 4. Comparison of our calculation (open 
symbols) with CHORUS data [TTl (full symbols) 
for neutrino (circles) and antineutrino (trian- 
gles) differential cross-sections on ^"^Pb at i? = 
35 GeV (upper plot) and E ^ 90 GeV (lower 
plot). 



Figure 5. Comparison of our calculation (open 
symbols) with NuTeV data ^21 (full symbols) for 
neutrino (circles) and antineutrino (triangles) dif- 
ferential cross-sections on ^^Fe at i? = 85 GeV 
(upper plot) and ii^ = 170 GeV (lower plot). 



